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Weak scale supersymmetry is a highly motivated extension of the Standard Model that has 
a strong degree of support from data. It provides several viable dark matter candidates: 
the lightest neutralino (a WIMP), the gravitino, and the axion/axino supermultiplet. The 
LHC turn-on is imminent. The discovery of supersymmetry at the LHC will go a long way 
towards establishing the nature of dark matter. I present arguments why mainly axion 
cold dark matter is a better fit for supersymmetric models than neutralinos. I also argue 
that Yukawa-unified SUSY GUT theories based on 50(10) with mixed axion/axino cold 
dark matter are extremely compelling, and present distinctive signatures for gluino pair 
production at the LHC. 

1 Introduction 

Astrophysical evidence for the existence of dark matter is now overwhelming, and comes from 
disparate sources: galactic clustering, galactic rotation curves, anisotropics in the CMB, mi- 
crolensing, large scale structure, to name a few. The dark matter clusters on large scales, and 
helps seed structure formation in the universe. While the identity of the dark matter particle, or 
particles, is unknown, we do know several of its properties: it must be massive, non-relativistic 
(cold or warm), electric and color neutral, and stable at least on cosmic time-scales. 

Of all the fundamental particles in the Standard Model (SM) of particle physics, only the 
neutrinos come close to having these properties. However, neutrinos are exceedingly light and 
engage in weak interactions: relic neutrinos would move at highly rclativistic velocities, and so 
couldn't clump enough to seed structure formation. Their measured abundance from WMAP 
analyses is only a tiny fraction of the universe's energy budget. Thus, the existence of dark 
matter is also evidence for physics Beyond the Standard Model (BSM). 

While there exist a plethora of candidate DM particles from BSM theories (examples include 
black hole remnants, Q-balls, sterile neutrinos, axions, KK gravitons, gravitinos, neutralinos, 
KK photons, branons and the lightest T-parity odd particle of Little Higgs theories. Two 
of these stand out in that they arise naturally due to very elegant solutions to long-standing 
problems in particle physics. These include the axion, which arises from the Peccei-Quinn 
(PQ) solution to the strong CP problem of QCDpQ, and the lightest supersymmetric particle, 
or LSP, of -R-parity conserving supersymmetric (SUSY) theories [3J. The SUSY theories solve the 
problematic quadratic divergences associated with scalar fields by introducing a new symmetry 
which relates bosons to fermions, thus giving scalar fields the milder divergence structure which 
is held by chiral fermions and gauge fields. SUSY also provides a means to unification with 
gravity, is an essential part of superstring theory, and in concert with Grand Unified Theories 
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(GUTs), receives some experimental support in the unification of gauge couplings under MSSM 
rcnormalization group evolution. 

In fact, the PQ strong CP solution and supersymmetry are in many ways made for each 
other, so these two schemes are not mutually exclusive, and both may well be right. In that case, 
the axion occurs as part of an axion supermultiplet, which contains along with the axion an 
i?-parity odd axino a, which may serve as LSP. In SUSY theories, the neutralino, the gravitino 
and the axino are all possible LSP candidates. In this talk, I will restrict my comments to 
supcrsymmetric dark matter (which includes axions and axinos), and comment on how it relates 
to LHC physics. 

2 SUSY WIMP (neutralino) cold dark matter 

In SUSY theories with neutralino CDM, the Xi is considered a natural WIMP candidate for 
dark matter. The neutralino relic density can be calculated by solving the Boltzmann equation 
as formulated for a FRW universe. Central to the calculation is computation of the thermally 
averaged neutralino annihilation (and co-annihilation) cross section. The fact that the relic 
density comes out approximately in the right ball-park is often referred to as the "WIMP 
miracle" . 

The paradigm model for SUSY phenomenology is called minimal supergravity (mSUGRA) 
or CMSSM. It is defined by just a few parameters: mo, miw, Aq, tan/3 and signal). In the 
mSUGRA model, the WIMP miracle is actually no miracle at all! The relic density turns out 
to be much too large over most parameter space. In fact, only in regions where the neutralino 
annihilation rate is highly enhanced will the relic density match the measured value. These 
regions are termed: the stau co-annihilation region at low mo, the HB/FP region at large mo 
where \i becomes small and we get mixed bino-higgsino CDM, the A-resonance region at large 
tan P where XiX? annihilation through the A- resonance is enhanced, and the (largely excluded) 
bulk region, where annihilation is through t-channel exchange of light sleptons. 

Direct production of WIMP dark matter at LHC (e.g. pp — > XiXi^Q is usually not inter- 
esting since there is no hard energy deposition for detectors to trigger upon. However, if SUSY 
exists, then LHC may be able to produce many or all of the other SUSY particles. The SUSY 
particle's subsequent cascade decays[3] should lead to collider events with high pt jets, high pt 
isolated es and /xs. Since each sparticle cascade decay terminates at the LSP (the putative DM 
particle), the SUSY events should also contain large missing transverse energy (E™ lss ) due to 
non-detection of the DM particles. 

At each point in mSUGRA (or any other SUSY model) parameter space, we can simulate 
LHC production of the entire array of superparticlcs, along with their cascade decays[3]. By 
looking for signals with high pt jets, isolated leptons and Evj? lss - beyond levels expected in 
the SM- we can test if a signal can be seen for an assumed value of integrated luminosity. In 
Fig. 11(51 [6], such a calculation has been made assuming 100 fb" 1 of integrated luminosity. LHC 
should be able to see the parameter space below the contour marked "LHC" , which corresponds 
to uig ~ 3 TcV when m g ~ rriq, or m g ~ 1.8 TeV when rriq ^> m g . We also show contours 
of direct WIMP detection rates and indirect WIMP detection rates via high energy neutrino 
detection at IceCube, or via detection of 7s, e + s or ps arising from neutralino annihilation in 
the galactic halo. 

It is noteworthy than in the DM- favored HB/FP region at large mo, the LHC can only cover 
a portion of allowed parameter space. However, in this region, direct detection via Xenon-100 
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Figure 1: Contours of direct and indirect detection rates along with reach of LHC and iLC for 
SUSY in the mSUGRA model for tan/3 = 55. 



or indirect detection via IceCube is likely. In the A-annihilation region- the large bump in the 
center of the plot- detection of halo annihilations via 7, e + and p is enhanced [7]. 

The enhancement of DD and IDD in the higgsino-like region is a general feature of a large 
assortment of models going beyond mSUGRA. In Fig. [2 we show predicted rates in models 
with a well-tempered neutralino[5]. The large cluster of models around 10 -8 pb shows that the 
next set of DD experiments can either discover or rule out an entire class of well-motivated 
SUSY models [5] . 



3 The gravitino problem for WIMP and gravitino dark 
matter 

A potential pit-fall in the mSUGRA model is known as the gravitino problem. If we assume a 
SUGRA-type model, with a TeV scale gravitino G, then gravitinos can be produced thermally 
in the early universe (even though they are never in thermal equilibrium). If the G is not 
the LSP, then it will decay into particle-sparticle pairs, and the sparticle cascade decays will 
contribute additional LSPs to the relic density. The relic density is too much if the re-heat 
temperature Tr ~ 10 10 GeV. Even if Tr is lower, the late-time gravitino decays inject high 
energy particles into the cosmic soup during or after BBN, which can destroy the successful 
BBN predictions which match so well with data. Detailed calculations [TU] show that one needs 

T R ~ 10 5 GeV (which conflicts with many baryogenesis mechanisms) or 10 5 < Tr < 10 9 GeV 
as long as ~ 5 TeV (the large ttIq suppresses the gravitino lifetime to less than 1 sec, so 
G —> f f decays occur at the onset or even before BBN starts). Since in SUGRA models the 
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Figure 2: Direct detection rates for SUSY models with a well-tempered neutralino. Each point 
represents a relic-density consistent model with f2^o/i 2 ~ 0.11. 



gravitino mass sets the scale for all the SSB terms, then we would naively expect all the SUSY 
particles to be at masses > 5 TeV (beyond LHC reach). 

One way out might be to assume the gravitino is the LSP. But this results in the gravitino 
problem in reverse. Then neutralinos or other sparticlcs in the early universe would decay 
with long lifetimes in particlc-gravitino pairs, and again disrupt BBN. Detailed calculations [TD] 



show that tuq should be 



1 GeV for m< 



100 



1000 GeV. Since rag sets the scale for the 



other sparticle, we would expect all the sparticles to have mass 
experimental limits. 



1 GeV, in contradiction to 



4 Mainly axion CDM in minimal supergravity model 

Another possibility is to assume some form of PQ solution to the strong CP problem. In 
the SUGRA context, we will add an axion supcrmultiplct to the model, which also contains 
an i?-odd spin ^ axino q|ll[ 112] . Then, axions will be produced as usual via vacuum mis- 
alignment [13] , and can contribute to the relic density. Their relic abundance depends on the 
PQ breaking scale f a , or alternatively on m a . A value of f a ~ 10 12 GeV, corresponding to 
m a ~ 10 -6 eV, would saturate the measured DM abundance. 

The value of the axino mass rria is very model-dependent: estimates range from the MeV to 
the multi-GeV scale[12]. In mSUGRA, if a is the LSP, then \i «7 can occur with a lifetime 
of order a fraction of a second: it is then BBN-safe. The a can also be produced thermally like 
the gravitino [14] . Its thermal abundance depends on f a , rria and Tr. Thus, DM would have 
three components: cold axions, cold axinos from thermal production and warm axinos from 



4 



Patras 2009 



neutralino decay. 

The scenario works well if ~ the MeV scale. Then, we find the mSUGRA model can 
have dominant axion CDM with a small mixture of warm and cold axinos[TS]. The resulting 
value of Tr is plotted in mSUGRA space in Fig. [3J The regions of mSUGRA space that are 
most neutralino dis-favored lead to the highest values of Tr > 10 6 GeV. This is enough to 
sustain non-thermal or Affieck-Dine leptogenesis! Thus, the most dis-favored neutralino DM 
regions are precisely the most favored axion/axino regions! Also, consequently, we expect quite 
different collider signatures at LHC in the case of mixed axion/axino CDM, as compared to 
neutralino DM. We remark here that calculating the amount of fine-tuning of the neutralino 
relic abundance also shows a preference for mixed axion/axino CDM over neutralino CDM[I^]. 




m„ (GeV) 



Figure 3: Contours of log 10 Tr (the re-heat temperature) in the mSUGRA model with mainly 
axion CDM, for tan f3 = 10. 



5 Yukawa-unified SUSY, mixed axion/axino CDM and 
the LHC 

SUSY GUT models based on the gauge group 5*0(10) are extremely compelling, since these 
models allow for matter unification into the 16 dimensional spinor of 50(10), give rise naturally 
to see-saw neutrinos, and yield automatic cancellation of triangle anomalies [17]. The simplest 
50(10) SUSY GUT models also contain t — b — r Yukawa coupling unification. Scans over SUSY 
model parameter space reveal that Yukawa unification only occurs for a very specific spectra: 
first/second generation scalars at the 10 TeV level, third gen. scalars and Higgs at the TeV 
scale, while gauginos are quite light, with ra§ ~ 300 — 500 GeV and m^o ~ 50 — 80 GeV [T5l[H)] . 
The neutralino relic abundance turns out to be much too high: fi^oft 2 ~ 10 3 — 10 4 , about 4-5 
orders of magnitude too much. The problem can be solved by invoking mixed axion/axino 
CDM[l9]. Then, Xi ~ * ^7 decays reduce the relic abundance by factors of 10 4 — 10 5 . The 
scenario works best if the CDM is mainly axions, with a small admixture of thermal and non- 



Patras 2009 



5 



thermal axinos [20]. Values of T R ~ 10 7 - 10 8 are possible, allowing for baryogenesis. And since 
rdsparticie ~ 10 TeV, we also expect ~ 10 TeV, thus solving the gravitino problem! 

In this model, the light gluinos lead to robust signatures at both the Tevatron[2l] and LHC 
colliders [52] , and the whole scenario should be largely tested within year 1 of LHC operation, 
or even sooner if a Tevatron analysis is performed. 
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